Wavefunctional approach to the bilayer v = 1 system and a possibility for a double 

non-chiral pseudospin liquid 
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We systematically discuss candidate wave functions for the ground state of the bilayer v = 1 as the 
distance between the layers is varied. Those that describe increased intralayer correlations at finite 
distance show a departure from the superflid description for smaller distances. They may support 
finite energy meron excitations and a dissipative collective mode in the place of the Goldstone mode 
of the ordered phase i.e. describe a vortex metal phase, or imply even an incompressible, pseudospin 
liquid, behavior. Therefore they describe possible outcomes of quantum disordering at finite distance 
between the layers. The vortex metal phase may show up in experiments in the presence of disorder 
at lower temperatures and explain the observed "imperfect superfluidity" , and the pseudospin liquid 
phase may be the cause of the thermally activated (gapped) behavior of the longitudinal and Hall 
resistances at higher temperatures in counterflow experiments. 

PACS numbers: 73.43.Cd, 73.21.Ac, 73.43.Nq 



I. INTRODUCTION 

The bilayer v = 1 quantum Hall (QH) system con- 
sists of two layers of 2D electron gases, each with a filling 
factor 1/2, that are brought together at distance com- 
parable to the average distance between electrons inside 
layers and the tunneling is negligible. The physics of the 
bilayer has been the focus of much of experimental and 
theoretical work. We will mention only two major exper- 
imental findings. First in the experiment of Spielman et 
al. Q there was a very pronounced ( "spectacular" ) zero- 
bias peak in the tunneling conductance between layers. 
Second in the experiments of Kellogg et al. @ and Tu- 
tuc et al. Q, in the counterflow setting of the bilayer 
(where in each layer the current is of opposite sign than 
the other) the resulting longitudinal and Hall resistances 
were dropping to zero in the limit of zero temperature. 
Both types of experiments signal a superfluid behavior in 
the pseudospin (relative between layers) channel of the 
system. An elaborated theoretical work of Moon et al. 
|4| described this behavior in the framework of an easy- 
plane ferromagnet (where the layer index represents the 
electron pseudospin degree of freedom) and ensuing XY 
model description. Therefore the superfluid behavior was 
expected. Usually it is described as a consequence of an 
excitonic condensate inside the bilayer where each exci- 
ton stems from the electron coupling to its correlation 
hole just opposite in the other layer. Nevertheless some 
discrepancies were noticed. In the first experiment of 
Spielman et al. the peak due to some unknown source of 
dissipation, is not as high and narrow as in the analog, 
Josephson tunneling experiments in superconductors. In 
the second experiment of Kellogg et al. and Tutuc et 
al. there were no Kosterlitz-Thouless transition signa- 
tures in the measured counterflow resistance. It seems 
merons (vortices) are liberated (due to some unknown 
cause) from their confinement in the expected superfluid 
down to very low temperatures. 

These findings point out that a careful investigation of 



quantum fluctuations or quantum disordering at finite d 
(distance between the layers) is necessary. One way to 
approach this question is to adopt the Laughlin approach 
to the fractional QH effect and look for approximate but 
very close to the ground states wave functions at finite 
d. That this is possible it was numerically demonstrated 
in Ref. @ for not large d. The functions that approx- 
imate the true ground states incorporate the effects of 
quantum disordering by allowing the presence of com- 
posite fermions (CFs) next to composite bosons (CBs) 
, another transformed electrons, characterizing the bose 
condensate of the assumed superfluid. CFs connect in a 
special way to the condensate to ensure and maintain the 
rigidity of the bose condensate. 

In this work we will systematically discuss candidates 
that we expect would approximate very well the ground 
state wave functions at finite d. Some of them we expect 
would be in a competition with the ones that are fer- 
romagnetically ordered and possess a Goldstone mode. 
Among them we will point out to one that describes an 
incompressible (in all channels) state. If the state were 
a true ground state, it would be a ground state of topo- 
logical phase that supports quasiparticles - merons with 
finite excitation energy. In fact this would be a realiza- 
tion of what is usually called a spin liquid phase with 
two kinds of semionic quasiparticles - i.e. a double spin 
liquid, with the only difference that here we work with 
pseudospin instead of spin. We will use interchangeably 
words double pseudospin liquid and pseudospin liquid, 
for short, denoting the same phase. This gapped phase 
might be the one that appears in counterflow experiments 
causing activated (gapped) behavior of the longitudinal 
and Hall resistances for a range of higher temperatures. 
Another competing possibility with possibly finite excita- 
tion energy for merons is a compressible version in which 
quantum disordering has allowed weakly coupled meron- 
antimcron pairs i.e. a vortex metal phase. This state 
supports a dissipative collective mode that comes in the 
place of the Goldstone mode of the ordered phase. The 
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physics of this state should be relevant for the expla- 
nation of the dissipation effects in the experiments, in 
which disorder, at lower temperatures, would dissociate 
the expected closed loops of meron-antimeron pairs in 
the topological phase. 

The paper is organized as follows: Sections II, III, and 
IV contain results of our paper, and Sections V and VI 
discussion and conclusions. Section II is an introduction 
to the wavefunctional approach to the bilayer, Section 
III in two different approaches with compatible conclu- 
sions describes the physics of a candidate wave function 
- vortex metal state, and Section IV introduces its mod- 
ification due to CF pairing - a pseudospin liquid state. 

II. WAVEFUNCTIONAL APPROACH TO THE 
BILAYER v = 1 SYSTEM 

If the distance between the layers is of the order or 
less the magnetic length - average distance between elec- 
trons inside the layers, interlayer Coulomb interaction 
will force the system, at total filling factor one and in 
the conditions of QH effect, to form the (111) state. The 
state is a simple generalization of the single-layer filling- 
factor-one case i.e. completely filled lowest Landau level 
(LLL) of one species electrons, precisely it is 

^Ul(zh Z l) = Y[( Z i] - Z J^]J( Z kl - Z ll)Y[( Z p1 - Z qi) 
i<j k<l p,q 

(1) 

where Zjf and Zi^ are two-dimensional complex coordi- 
nates of electrons in upper and lower layer respectively 
and we omitted the Gaussian factors. If layers are far 
apart each one will be a separate system at filling fac- 
tor one-half. Numerous studies show that the ground 
state at that filling factor is a generalization of the Laugh- 
lin construction which includes a single Slater determi- 
nant of noninteracting particles in zero magnetic field, 
i.e. 

*i/aM =T{<S> f ( W ,w)Y[(w t T -™ jT ) 2 } (2) 

i<j 

where $/ is the determinant and V represents projec- 
tion to LLL. In Refs. |6|, |7| it was shown that relevant, 
underlying composite particles in this case are CFs. On 
the other hand the relevant, weakly interacting compos- 
ite particles in the (111) case are CBs H,0). 

Now let us start from (111) case, increase the dis- 
tance and introduce one- half correlations in the (111) 
state in a minimal way, preserving (111) intercorrelations 
of newly introduced CFs with all other remaining CBs. 
This means that though we are perturbing (111) state, 
we are assuming its inherent rigidity. The wave function 
that describes this is 

m bbf =VA { Y[(zii - z n )Y[(z kl - z n )Y[(z p] - z ql ) 

i<j k<l p,q 



x $/(w T ,u; T ) JJ(w lT - Wj] f 

i<j 

x ®f{u>i,wi) Y[(w kl -w n ) 2 

k<l 

i,j k,l 

x Y[( z t l~w q ^Y[(z ml -w nl )} (3) 

p,q m,n 

where we omitted the Gaussian factors, V denotes the 
projection to the LLL, A - antisymmctrization between 
bose and fermi variables in each layer separately, and <f>/s 
are Slater determinants of free waves. In the thermody- 
namic limit the flux - number of particle relationship is 

iV* = N b] + N H + JV> t + Nfi 

= 2N n +N b1 +N bL (4) 
= 2N fl +N b1 +N bl . 

Consequently we must have Nff = Nf± but there is no 
constraint on the relative number of bosons. 

In fact, once we adopt the hypothesis that with in- 
creasing d (distance) CFs are slowly nucleating we are 
left, because of the flux-counting arguments, with only 
two (simple, ala Laughlin) possibilities for interpolating 
ground state. The second possibility has instead of the 
last two lines in © for the intercorrleations between the 
bose and fermi parts, the following expression, 

IJ(z iT - w n ) 2 Y[(z kl -wn) 2 . (5) 

i,j k.l 

These intercorrelations are of the fermi part type, i.e. 
they are intralayer correlations, a consequence of possibly 
more important intralayer Coulomb interactions. 
In this case we have 

N<& = 2N f] + 2N b] 

= 2N n + 2N H 

= 2N n +N b1 +N bL (6) 

- 2N n + iV bT + N bl 

i.e. both bose and fermi numbers are constrained to 
7V /T = N fl and iV 6T = N bl . 

In the small particle (5+5) numerical study in [fj the 
first possibility had much larger overlaps, presented in 
0, with the true ground state as it was evolving with 
distance (and the ratio between bosons and fermions was 
changing) than the second possibility (Eq.©). The over- 
laps of the first possibility were slowly decreasing with 
distance so that no definite conclusions can be drawn 
about the precise evolution of the system near the tran- 
sition region. But for smaller d the relevance of the first 
possibility for the evolution was established. 

We will not discuss the states that we naturally ex- 
pect would interpolate between these two limits, small d 
(Eq.©) and near the transition (Eq.JSJ) limit. In these 
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states some of the fermions would connect via bose (111) 
type intercorrelations to the bose part and some via fermi 
type. 

It is no wonder that the first possibility is more relevant 
for smaller d. It also allows the imbalanced (Ni ^ N2) 
situation, which is required due to the theoretical and 
experimental insight gained 0, 0, HH ITU ] about the exis- 
tence of the Goldstone mode connected with the particle 
number difference. (The projected to definite particle 
number picture, which we here discuss, must incorporate 
that physics.) Later we will show that the state is incom- 
pressible (in the charge channel), at least in the scope of 
a suggested Chern-Simons picture. 



A. Chern-Simons approach 

In this section we will employ the classical CS ap- 
proach, described in |8} and |7j with its usual RPA (ran- 
dom phase approximation) to find out the linear response 
of the system that supports t he g round state in Eq.®. 
This approach is advanced in [l2j but at the same time 
there, a conclusion is drawn that the classical CS ap- 
proach in the RPA allows us with not much work to find 
out qualitatively the physics of the linear response. Here 
we have still another problem; there is no obvious way 
to implement antisymmetrization procedure of the wave 
function in Eq.©. The antisymmetry of electronic wave 
functions comes naturally and automatically only in only 
CB or only CF theories. Therefore to simplify matters 
from the start we choose the classical CS approach and 
neglect the antisymmetrization requirement. 

Without the antisymmetrization the wave function in 
© would be the ground state (in the RPA and when 
the projection to LLL is neglected) of the following CS 
theory 

C = 5^{»*t(fl,+taS-*4>-«rB )*a 
+ ^-Y^^Adk + ial-iA k -iaB k )H a } 

k 

+ ^2{i<f>l(d + iao - iAo - iaBo)<f>„ 

a 

/m z — ' 



1 1 



^ 2vr 2 V ' 



\Y,J d?5 Pa (r)V a {r-7)5p a {?) (7) 



dr'5p ] {r)V e (r-r')8p i (r' 



where a is the layer index, taking f and { values as a 
variable, ^ a and are CF and CB a-layer fields re- 
spectively, and V a and V e are intra and inter-Coulomb 
interactions. <5p CT 's are given as sums of CBs and CFs, 



i.e. Sp a — Sp^ + 5p^ , as the most natural choice for the 
electron density in this distinguishable picture. External 
fields, A and B, couple to charge and pseudospin (up 
minus down) degrees of freedom respectively. The CS, 
gauge fields, a-p and a^, are related to a as 

_ a T + a l 



and in this way, in the mean field approximation, repro- 
duce the relations encoded in the ground state given by 
Eq.@. 

The antisymmetrization can be implemented by the 
following constraint, 



S-S 



N* ( N* 
2 [ 2 



1) 



(8) 



for each layer separately so that No- denotes the number 
of electrons in the layer. S denotes a generalized spin of 
the layer obtained by integrating over the volume of the 
system of the following field density, 

$(r) denotes a spinor for which we have, 



where Ub and Uf a are the CS unitary transformations, 



and 



+Pfa(r') + Pf-a(r')) 



UfAr) = 

e XP U^ar 9i r-r' ){ 2p M r' ) + 



(9) 



Pba{r') + Pb- a {r')) \, 



(10) 



$ a and ^cr, are already introduced bose and fermi field 
respectively, and a arc the usual Pauli matrices. The idea 
behind the constraint in Eq.© is simple; it uses the fact 
that spin i particles, N a of them, must necessarily, in the 
states for which the constraint in Eq.JSJ) is true i.e. have 
S-S at its largest possible value, be completely symmetric 
in the spin space and necessarily antisymmetric in the 
real space. 

We still have to fix S z , which is the second constraint, 
i.e. the number difference between CBs and CFs in the 
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layer, to project the artificially introduced spin \ problem 
to the reduced Hilbert space and our problem. 

What we can immediately notice is that the con- 
straint in Eq. (jSJ) will introduce the terms that interchange 
bosons and fermions at different points instantaneously 
which conforms to our idea of indistinguishability. Al- 
though it is easy to formulate it is very hard to im- 
plement the constraint. Only, maybe, if we have over- 
whelming number of bosons or fermions we will be al- 
lowed to neglect the constraint as it is usual in the case of 
bose fluids with the natural decomposition into classical- 
macroscopic, condensate part and normal part. 

Nevertheless there is a much deeper reason that allows 
us to neglect the antisymmetrization requirement. The 
reason is that, just like in a hierarchical construction and 
as it will be much more clear later, CFs represent meron 
excitations (meron-antimeron pairs, see Appendix) that 
quantum disorder the (111) state. As it is usual when 
we discuss the dual picture of the FQHE [l3( , we do not 
extend the antisymmetrization requirement to the quasi- 
particle part of the electron fluid. 

We want to prove that indeed the static density- 
density response that follows from the Lagrangian in 
Eq. (O shows the incompressibility of the system with 
the ground state in Eq.© in which the antisymmetriza- 
tion is neglected. The charge of the system is related to 
5a as 



5p = 5pi + 5p i = 



ikSa 

~2tT 



static (uj = 0) case we are left with the following charge 
part, decoupled from the pseudospin part in the RPA, 

C e ff = Koo(Sao) 2 + Kn(Sa) 2 + —a ika 



2tt 

2m 1 "" 1 2 (2tt)2 



V c (k). 



The integration over 5p b c gives the constraint ao = Ao 
and consequently the integration over 5a gives for the 
(static) density-density response, 



2p, fen^fe) > 

where K\\ = — 12 ^ m Q. As long as there is a finite den- 
sity of bosons (p fc ) the system is quantum Hall, incom- 
pressible (and the above expression vanishes in k — *■ 
limit). Even the presence of composite fermions i.e. par- 
ticles which we would naively consider to represent gap- 
less degrees of freedom in the system, does not lead to 
compressibility. A small fluctuation in the charge den- 
sity (5pf = 5p^ + 5 pi ) of the composite fermions lead to 
a fluctuation in the bosonic charge density for which we 
know we need a finite amount of energy to create § . 



where, as we work in the transverse gauge, \5a\ = 5a is 
the transverse space component of the vector. For the 
response calculations (with Bo = B = 0) we will adopt 
the conventions introduced in 0. Integrating out the 
quadratic terms in fermionic fields to the second order 
in gauge fields (RPA) Q , then introducing density-angle 
variable for the bosonic fields and expanding the action 
in them neglecting the amplitude-density derivatives § , 
we get 

^ff = E^ X oo(K) 2 + ^n(^) 2 ] 

a 

+ Yl h Ba (-k){iu®a{k) - 5aZ{k)) 

a 

2 (2ir) z 2 2ir 

where Ve (k) = MM+M Vs (k) = WW , K Q0 = 

Kqq = Kqq and K\\ = k} x = K xl are density-density 
and current-current response functions of a free fermionic 
sustem respectively @, ©f, ©1 are bosonic angle vari- 
ables and p b averaged bosonic density per layer. In the 



III. A STATE WITH POSSIBLY DECONFINED 
MERONS 

The second possibility, state 
^bff = VA{ Ui<j( z H - z jl)Tlk<l( z H ~ z n) 

n^^T-^T^nM^i-^u) 2 }^!) 

should be relevant for the transition into two decoupled 
Fermi seas region. In the following section we will explore 
its properties; in the first part we will show that in this 
state possibly a finite energy is needed to excite meron, 
and in the second part we will show that the state is 
compressible in the pseudospin channel and nearly sup- 
ports a gapless pseudospin mode. The first property tells 
us that in this state merons may be deconfined relative 
to the (111) and Eq.© states in which merons as vor- 
tices of the pseudospin superflud are confined. If it were 
not for the second property, the state would possibly de- 
scribe a quantum Hall, topological phase with four kinds 
of gapped, meron quasiparticle excitations. 



A. The screening of meron 

In the following we will show that the screening charges 
of the meron excitation in the plasma analogy of the state 



5 



in Ea. l|llf> are localized, not of long range, and may not 
lead to the usual logarithmic divergence with the size of 
the system of the energy required to excite a meron. 

Namely we will use an effective expression from 
Ref.^H for the meron excitation that in the (relative 
change in) density calculations for the meron state leads 
to results valid in the longwavelength limit. The effec- 
tive expression is the one reduced from the well known 
expression in the (111) state Q, and in the Ref. |3| it 
was shown that leads to the logarithmic divergence in the 
energy to excite meron in both (111) and the first mixed 
state (Eq.©) that describe the pseudospin condensate 
at small d. Therefore the effective expression gives the 
expected behavior of a meron excitation in these states. 
As in the Laughlin quasihole construction, the meron ef- 
fective construction multiplies the ground state, \& G . Ex- 
plicitly the construction is 



n 



Rt 



w 



C 



2|z. 



it 



w 



(12) 



where C is a constant, found to be equal to C = 0.80 in 
the case of the (111) state. We emphasize that Eq. <|12[) 
is an effective expression, modeled for the task of density 
calculations, and in no way an expression that would be 
valid in the short distance limit or stand for a meron 
construction in the lowest Landau level. 

In the plasma language of the Ref.0], and as an 
interpretation of the squared norm of the construction 
(Eq.(^2J0 we have that an impurity at point w connects 
via the interaction C/\zii \ to the f particles of the plasma 
defined by the state, ^> . To obtain the charge contribu- 
tions (| and I) far away from the center of excitation 
in the second state (Ea.^lip). i.e. if we take in Eq. l|ll|) 
that ^> is the second state, we use the sam e ty pe of 
approximations explained in Refs. 0] and [lif. We 
assume that the classical partition function of the clas- 
sical system that is defined by the square norm of the 
second state has the property of screening, and there- 
fore that the essential physics including the description 
of the screening charges of impurities can be found by 
summing so-called chain diagrams. The contributions to 
the screening charges of an impurity can be easily vi- 
sualized as a type of chain diagrams shown in Figs. 2 
and 3 which connect impurity (w) to the probing point 
(r) on the right-hand side. In the Figures the straight 
line denotes the - interaction, i.e. - in the momentum 

r ' q 

space by which the meron in Eq. p2|) connects to | parti- 
cles, the wriggly line denotes the ln{r}(~ ^) interaction 
i.e. the Coulomb interaction of 2D plasma, and there 
are two kinds of vertices. For bosonic quasiparticles we 
have a vertex that equals to their densities, n-f or nj., 
but for fermionic quasiparticles the value of the vertex is 
the static structure factor of free Fermi gas i.e. S|(g) or 
si(q) for which we have s^(q) = s^(q) ~ q in the long- 
wavelength limit. To calculate the screening charges let 
us introduce two infinite sums; the one for j particles is 



FIG. 1: The infinite sum of Ea. fHty 



FIG. 2: The up density screening 



depicted in Fig. 1 and equals 
4s T U 2 n? 



1 - 2sfV 



-2Vn 2 , 



(13) 



where V ~ Jy, and similarly we have JJ. for j particles. 

Using these infinite sums we are able to write the effec- 
tive interaction between n-\ vertices in a compact form, 



V(n T ,n T ) 



i Vr. 



l-V(n+K+]j.) t 



l-V(n+fr+^) 1-2-H- 



•T 



(14) 



In Eq. (|14(l we also stated the value of the interaction in 
the longwavelength limit. The introduction of the infinite 
sums comes naturally because of the type of intercorre- 
lations that exist in the electronic wave function. In this 
way we can express the effective interaction between a 
bosonic and fermionic f vertex, V(n.-]-, s-f), as 



V(n T ,s T ) 



-,,/ \ 2s r v q^o _ 



(15) 



with as the value in the longwavelength limit. Similarly 
we have for the effective interaction between fermionic f 



FIG. 3: The down density screening 
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vertices, 



v(*mt) 



( St ) 2 2V + 
V(sT. n J i- 2sT v — ► -s T (t?). 



(16) 



Now we can combine all these effective interaction expres- 
sions to find out the screening charges with up pseudospin 
of the meron excitation in Eq. (|12|) . The first contribution 
(pbb in Fig. 2) is due to the meron connection to | bose 
quasiparticles and, if we denote the direct interaction to 
t bose quasiparticles by V„, 



it is 



(17) 



Then the second contribution (pbf in Fig. 2) is through 
the direct interaction to up bose quasiparticles that ends 
up with up fermionic quasiparticles; it is 



y m V("T.«i)«z3 , 



(18) 



The third contribution (pfb in Fig. 2) comes from the di- 
rect connection to up fermionic quasiparticles that ends 
up with up bose quasiparticles. This contribution is sim- 
ilar to the previous one and has the same limit. The last, 
fourth contribution (pff in Fig. 2) connects through up 
fermi quasiparticles again up fermi quasiparticles, and 
equals 



h s l) 



*1 



(19) 



A comment is in order here. Although our plasma ap- 
proach is straightforward and leads clearly to the local- 
ized meron screening charges it does not automatically 
lead to an overall finite meron excitation energy because 
we can draw a conclusion only for the charging part of 
the energy (charge (f and j) difference squared). But it 
certainly signals a possibility for meron deconfinement. 



B. The pseudospin degrees of freedom 

Now, in the scope of already introduced Chern-Simons 
theory, we want to investigate the pseudospin channel 
of the state in Ea. Hll(l . In particular we want to know 
wheather the pseudospin degrees of freedom are com- 
pressible and support a gapless mode. For this we need a 
variant of the Chern-Simons theory introduced in Eq. J7|) 
for the first state introduced in Eq.©. In this case we 
have four gauge fields af a ,a Bcr ; i = 0, 1, a =f, J, 



V x a Fa = 27r(2p Fa + 2p Ba ) 
Vxa Ba = 2tt(2p Fct 



p B ^+p Bl 



(22) 



acting on ^ a and fermi and bose fields respec- 
tively. In fact, as can be easily seen, we have only 

BT_j_ Bl FT I PI 

three independent gauge fields a c — - — ±2 — — - — ±2 — 



and 



a ' —a 



and a. 



2 2 

and, therefore, 



fs 



ika\ 
2tt 



= P 



(23) 



This contribution is again equal to zero according to 
Eci. (|16fl in the longwavelength limit. 

For the meron screening by down quasiparticles we 
need V(n|,n;) interaction, and similarly to the same 
pseudospin interaction we get 



1-V(n+1H4) 1-2-&. ' 



(20) 



where (8 is f|" or JJ. and n a is n-f or n\ whether we attach 
(g> to the right or left side of a diagram respectively. In 
this way we get, in the longwavelenght limit, 



V(n T ,nJ?=5o. 



(21) 



Because of this any contribution of the down screening 
charges in which participate bose quasiparticles is equal 
to zero. In fact V(n-f, n{) participate in each of four pos- 
sible contributions, see Fig. 3, leading us to the con- 
clusion that both, p-f and screening charges tend to 
zero in the longwavelength limit. Therefore the screen- 
ing charges are short ranged and localized just like in 
the case of any one-component quantum Hall system - a 
topological phase. 



pK 



(24) 



Then, similarly as before, we can write down the effective 
Lagrangian for the pseudospin part, 



C ps = ±a° ika{ s + ±a f s zk(af - a{ s ) - \(±) 2 V s \a{ 



oo -t 
lm <\5p b s \ 2 



4 ft? 



2V27T. 

al\ 2 Kn + Sp s B a 
4 s 6p b s -%\a{r 



(25) 



We first integrate out Oq S which gives us the expected 
constraint on the pseudospin density of bosons, Spf . Us- 
ing this constraint and then integrating out a{ s and af 
we get for the pseudospin density - density response the 
following expression, 



C eff 



1/ i: \2n2 



W* + (±) 2 V s -2K n -±=^ 



(26) 



where 



W 



f 1 rncu 2 \ 



( — ) 2 - 

* 2?r ' 



2K„, 



1 mu 2 , 2p^ 

2 p(,(2ir) 2 ' m 



(27) 
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First in the limit u> — > (and then k — > 0) we see that 
the system is compressible. Second, taking into account 
that 



k 2 2p f co 

h l — - — , 

12nm kp k 



(28) 



we see that the pseudospin gapless mode, in the case of 
the (111) state simply lj° = \J 2 ^ L ~p b k, does not exist as 

an eigenmode in this case. It is nearly so if we have in 
mind that the fraction of bosons in this state is to be 
considered small, and so is the imaginary part of K\\ 
when for the frequency we consider the one that takes to 
zero the real part of the denominator in C e ' ' . 

In fact also for the state in Eq. © we find that there is 
a pseudospin eigenmode with a dispersion relation lu = 
oj°(k)+ick 3 where c is a constant. Therefore it is slightly 
dissipative what we do not expect from the Goldstonc 
mode |4|. In the following section we will consider the 
variational constructions, Eq.(J3J) and Ea. (|ll|) . with p - 
pairing of composite fermions, introduced in Ref. |l6( 
, that cures the dissipation problem of the construction 
in Eq.lO and leads, as we will show, to incompressible 
behavior also in the pseudospin channel of the state in 
Eq.jnj. 



IV. COMPOSITE FERMION PAIRING AND A 
POSSIBILITY FOR A PSEUDOSPIN LIQUID 

In Ref. p-wave CF pairing was proposed as a way 
of lowering ground state energy of the state in Eq.©. 
That was explicitly shown on a basis of Monte-Carlo cal- 
culations in which the paired states are excellent varia- 
tional ansatzes with respect to the true ground states. In 
the context of the phenomenological Chern-Simons the- 
ories based on the proposed wave functions it is not hard 
to show that in the case of p-wave CFs pairing the state in 
Eq. acquires a pseudospin mode without an imaginary 
term which must be then a Goldstone modepredicted by 
the theory of the ordered state for small d Q • 

The calculation begins by noting that p-wave pairing 
of CFs is simply a condensation into " 11-1" state by the 
way of the Cauchy identity, 



dct{ 



Wpf—Wql 



(29) 



In fact, as stated in Ref.0| , the p-pairing they found 
is with the pairing function g(w) ~ i.e. a p x — ip y 
instead of p x + ip y pairing, and in the context of the 
Chern-Simons theory that means that the gauge fields 
on CFs acquire negative sign, so that in the end we have 



Zfc Fa _ ife Ba — F I F 

2n ai ~ 2ir ai -P° +p -° 



Pa + P B o 



(30) 



Then it is easy to show, similarly as before, by integrating 
out dpf and Spf , in the pseudospin channel, that the 



Goldstone mode has dispersion ui(k) — y (Ph + Pf)'k- 
We get the same result if we consider p x + ip y pairing 
instead of p x — ip y . 

On the other hand if we consider the state in Eq. Hll|) 
with fermi type intercorrelations between the fermi and 
bose part and introduce the p x —ip y pairing between CFs, 
in the CS language we have 



ik 
2^ 



and 



ik 
2^ 



a?°=pZ+p B a + 2p- : 



(31) 



(32) 



i.e. § a( s = pf - pi and § of = pf - p? . Then the 
pseudospin part of the effective Lagrangian is 



reff = 



l™*\ Sp B l 2_ a fs Sp B_-£ la fs l2 



+ 



1 mu) 
4 p F k' 2 



\Sp 



Fl2 



«S*pf-|rlafl 



-\{±) 2 V s \a\ + 4 s ? ~ IK + 4 S )Bo- (33) 

In few steps, by reducing L e JJ into an effective, quadratic 
expression in B a we can find out the density-density cor- 
relator. It vanishes in the k — > limit and signals that 
the state in Ea. (|llf> in which CFs pair is an incompress- 
ible state. It is our expectation that in this state exist 
four kinds of merons, characteristic also to the ordered 
pseudospin state. Therefore the state, if a ground state 
of an electron system, should represent a quantum Hall 
i.e., in the low-energy limit, a topological phase. 

Even in existing numerical work we find a support for 
our expectation. In Ref. > i n t ne data that repre- 
sent the excitation spectrum of the bilayer v = 1 system 
at d = 1.51b Qb is the magnetic length), i.e. in the 
transition region, obtained by exact diagonalization on 
a torus, we can see signatures of the nearby topologi- 
cal phase. Namely, in the low-lying spectrum dissoci- 
ated from the Goldstone mode excitations, exists a four- 
fold degenerate energy level that we expect represents 
expected four ground states on the torus of the topo- 
logical theory [lg- (We do not interpret these states as 
ordered spiral states as in Ref. 17] because no nearby, 
low-lying excitations can be seen in the existing data.) 
Because of the importance of the intralayer correlations 
in the transition region, the mixed state with fermi type 
intercorrelations between the bose and fermi part should 
compete with the true ground state and is very relevant 
to the physics of the region. The topological theory in 
question is 1^(1) ® ^2(1) because it supports two kinds 
of semionic quasiparticles 01 • (Meron fractional statis- 
tics is semionic.) On the other hand the counterflow ex- 
periments 0j do not support a "perfect" superfluid 
scenario with a Kosterlitz-Thouless transition. Instead 
the data on the counterflow longitudinal resistance show 



8 



superfuid 
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two Fermi seas 






*~ 
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FIG. 4: The possible phase diagram of the bilayer as a func- 
tion of the distance between layers in the presence of disorder 
at zero temperature. 



Eq.© can be viewed as creation of two closely spaced 
paires of merons. Each pair consists of two merons of the 
same vorticity but opposite charge. On the other hand 
the creation of (2 + 2) CFs in the condensate of CBs as in 
the construction in Eq. leads to meron pairs that are 
not closely spaced inside each pair. (The construction in 
Eq. (|ll|) suits the increasing intralayer correlations with 
d.) This motivates a picture of the intermediate phase 
as either a condensate of closed meron-antimeron loops 
of any size - a pseudospin liquid [ly (section IV) or a 
soup of dissociated meron-antimeron pairs (loops) i.e. a 
vortex metal (section III) . The presence of disorder may 
stabilize the vortex metal phase as in Fig. 4, or maybe a 
phase separated version [ill ]. 



an activated behavior, for a range of higher temeratures, 
very similar to the usual data of a quantum Hall phase. 
Therefore it might be that due to an increased impor- 
tance of the intralayer correlations in the counterflow ex- 
periments, the topological phase stabilizes and with it a 
gapped behavior even in Hall resistance 2j, Q . Because 
of the non-chiral flow of the currents in the experiment 
the relevant topological theory should be [^(l) <8> ^2(1) 
double non-chiral pseudospin liquid. The theory is invari- 
ant under combined time reversal and Z2, an exchange 
of layer indecies, operations implied by the experimental 
setup. Furthermore it is well-known that the zero-bias 
peak in the tunneling conductance is not as of the usual 
Josephson effect in superconductors due to dissipation. 
The reason for this should also be found in the physics 
of the state in Eq. |(TJ that incorporates the effect of the 
increased intralayer correlations in the transition region 
in which the experiment occured. 

V. DISCUSSION: PHASE DIAGRAM FROM 
THE WAVEFUNCTIONAL APPROACH 

When we take into account what we found out in previ- 
ous sections, in the scope of the wavefunctional approach, 
the phase diagram, in the absence of disorder, may well 
have an intermediate phase between the superfluid phase 
and the two decoupled Fermi- liquid like phases. The 
wavefunctional approach tells us (sections III and IV) 
that in the intermediate phase the pseudospin stiffness 
may go to zero (meron deconfinement) but the density of 
composite bosons stays finite. The density of bosons dis- 
appears at the transition to the two decoupled CF Fermi 
seas. 

This scenario of the influence of quantum disordering 
in a superfluid is often described as a result of creation of 
vortex-antivortex pairs (loops) that cause phase fluctua- 
tions. The phase fluctuations may cause the disappear- 
ance of the superfluid long range order, but may leave 
the ordering amplitude (boson density) nonzero. In the 
Appendix we will argue that creation of (2 + 2) CFs in 
the (111) condensate of CBs as in the construction in 



VI. CONCLUSION 

By systematically investigating possible candidates for 
ground state wave functions for the bilayer v = 1 system 
we reached the conclusion that the ordered state that 
supports a Goldstone mode in the pseudospin channel 
for general d can be described as a mixture of CBs and 
p-paired CFs with dominant bose (111) type intercorre- 
lations between the two components. Because of the in- 
creasing intralayer correlations with d, a state that com- 
petes with the ordered one can be described as a mixture 
of CBs and p - paired CFs with dominant fermi type in- 
tercorrelations between two components. This state that 
describes a topological phase with four gapped meron 
quasiparticles may cause (already observed) activated 
behavior (decaying exponentially with a gap) behavior 
of the counterflow longitudinal and Hall resistances at 
higher temperatures. At lower temperatures, because of 
the presence of disorder, a compressible in the pseudospin 
channel version of this state, a vortex metal, may come 
out and cause the observed phenomenology of the bilayer 
"imperfect superfluid" . 

There is a vast literature on the bilayer. We are plan- 
ning a separate publication on the features of the vortex 
metal phase and comparisons to the previous, experi- 
mental and theoretical work. Here we will mention only 
theoretical work that like ours stresses the importance of 
the interactions and quantum disordering. In Rcf. 20] an 
excitonic square lattice solid phase was proposed. Here 
we discussed only translationally invariant, homogenous 
possibilities and inhomogenous states, Ref. 20], may also 
come as candidates. (The inhomogenous, phase sepa- 
rated version of the mixed CB-CF states was first pro- 
posed in Ref.[2l| before the homogenous kind In 
Ref. |22j] interactions only drive a phase transition into 
a state with algebraic (quasi-)long-range superfluid or- 
der. The numerics, without impurities and not biased by 
boundary conditions and special number of electrons to 
the Wigner crystal formation, Ref. 0] > tells us that only 
the true superfluid order exists in the pertinent experi- 
mental region. Therefore we want to stress again what 
our work suggests. Quantum disordered states ("dou- 
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ble pseudospin liquid") are nearby the true superfluid 
ground states for distance between the layers pertinent 
to the experiments. We need also impurities to stabi- 
lize a version of these disordered states, "vortex metal 
state" , at lower temperatures. Interactions caused quan- 
tum disordering is necessary but not sufficient condition 
for the explanation of the dissipative phenomenology of 
the experiments. 
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part of this work was completed. She also thanks M. 
Levin, N. Read, Z. Tesanovic, and X.-G. Wen for discus- 
sions. The work was supported by Grant No. 141035 of 
the Serbian Ministry of Science. 



APPENDIX A 

In this Appendix we will give arguments that the four 
(2 + 2) CF inclusion into the (111) state, i.e. 

S T {(zi T - z 2T )(exp{ifcz u } - exp{ikz 2l })} x 

s li( z il ~ z 2 i)(exp{ikz n } - exp{ikz 2l })} * m , 

(Al) 

where <Sf and Si are symmetrizers inside each layer, and 
we omitted for the sake of simplicity the overall projec- 
tion to the LLL, corresponds to two meron pairs, of op- 
posite (up and down) vorticity. Inside each pair merons 
have the same vorticity but opposite charge. Therefore 
the two meron constructions that can act on the (111) 
state and produce such a pair is (see Eq.^20 

= n.^^exp^^f^} x 

\z^-w 2 \ eX P{£.7 2\z,°-wo\ } ' 



■ij 2\z jr -w 2 \ 



(A2) 



and analogously for the down pair. These expressions are 
effective i.e. valid in the long distance approximation. 
Nevertheless in this Appendix we will take that they are 
qualitatively correct even for shorter distances. Then in 
the quasiparticle ("fractional statistics") representation 
[23| the wave function that describes the two merons in 
Ea. (|AT|l is 



^(wx,w 2 ) = (wx - w 2 ) 2 f(w 1 ,w 2 ), 



(A3) 



where due to the mutual semionic statistics between 
quasiparticles we have the difference, (wi — w 2 ), to the 
power — i, and f(w\,w 2 ) is a symmetric function of co- 
ordinates, in our case: 

f(wi,w 2 ) = (expjifcwi} — exp{ifcu?2})-7-!r ~rr- (A4) 

\wi-w 2 \ 

There are no Gaussian factors because, when calculated, 
the interaction of a meron with a positive background of 



the corresponding plasma, 
1 



d 2 z 



, ex P {-^}, 
z — w\ 2 



(A5) 



is a bounded function of w, and when expontiated gives 
a factor that weakly depends on w. 

To get the wave function in terms of electronic coordi- 
nates we have to calculate [23| 

J d*W 3 J d2 W4 _l_ /(w3;W4 )^ (w3) ^ (u;4) tf m . 

(A6) 

The combined exponentials in the two meron construc- 
tion can be expanded as in the following, 

c 



i i c; v 1 {zi—w-L)s 

2 {z i -iu 1 ) 2 \z i -w 1 \ 



r I o (zj-wi)S -| 

2 l^il (z z -w 1 ) 2 \z i -w 1 \ ~r ° (z t -w 1 ) 4 \z i -w 1 \ > 



2! V 2 ) (zi-W!) 6 



{zi — w\)8-{zj —w\ )<5 

M\S\ 3 ) 



2! V 2 > ^i^j (z l ~w 1 ) 2 \z i ~w 1 \(z J -w 1 ) 2 \z j -w 1 \ 



(A7) 



where 5 = w 2 — W\, and for the sake of clarity we su- 
pressed the up arrow. To the same order of accuracy we 
can rewrite the above expression as 



1 _L C V-> (zj-W 2 )S 

' 2 Z^i (zi-w 2 

-L(£\2y^ 

2! V 2 / L-ii (z i ~W2) 2 \zi~w 1 \(zi-w 1 ) 2 \zi-w 2 \ 



-w 2 ) 2 \zi—w 2 \ 
(zi-Wi)S-(zi — w 2 )5 



2! V 2 > t-ii^j (z 



(zj— wi)S-(zj— w 2 )8 



2 ) 2 | Wi - W 1 | ( Zj - W 1 ) 2 | Zj 



U 2 \ ' 



(A8) 



The contribution of the first two terms must be negligible 
or zero after the integration over wi and w 2 which brings 
averages of random phases with respect to electron dis- 
tributions encoded in \&m. The contribution of the third 
term, after picking sigularities at Zi, is identical to zero, 
but the contribution of the fourth term, after picking up 
singularities at z^ and Zj-\ is 



f(z lh z n ) ■ \zn - z n \ 

n (z kr -Z ir ) -i-r (zif—Zjf) 
k =£i \z k ^-z ir \ llljtj \z n -z jr \ ■ 



(A9) 



for the up part of the construction in Ea. (|A6(l . We get 
a similar contribution for the down construction, and 
therefore, up to some phase factors, we get Ea. l)Aljl . The 
phase factors are there to ensure the right flux count 
through the system, the problem we neglected by writing 
our wave functions (Eq.(J3J and Eg. Hill) ) in the thermo- 
dynamic limit. 

By examining (2 + 2) CF construction in which Fermi- 
type intercorrelations are dominant, i.e. the one as the 
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construction in Eg. i|ll|) for (2 + 2) CFs, we can find out higher order terms of the expansion in Ea. (|A7|l are rele- 
that the small distance approximation as in Ea. HA8|) is vant, and therefore two merons are well separated in this 
not applicable. Moreover its form suggests that much construction. 
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